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We look at various correlation functions, which include those that involve both the velocity and 
the vorticity fields, in two-dimensional (2D) isotropic homogeneous unforced turbulence. We adopt 
the more intuitive approach due to Kolmogorov (and subsequently. Landau in his text on fluid 
dynamics) and show that how the 2D turbulence's results, obtainable using other methods, may be 
established in a simpler way. Same method is used to calculate some third-order structure functions 
for quasi-geostrophic (QG) turbulence for the forward cascade of pseudo-potential enstrophy and the 
inverse energy cascade in quasi-geostrophic turbulence. These results motivate us to study the two- 
point third order structure function in the context of the two-dimensionalisation effect. Consequent 
studies enable us to give a reason for the inverse energy cascade in the two-dimensionalised rapidly 
rotating three dimensional (3D) incompressible turbulence. For such a system, literature shows a 
possibility of the exponent of wavenumber in the energy spectrum's relation to lie between -2 and 
-3. We argue the existence of a stricter range of -2 to -7/3 for the exponent in the case of rapidly 
rotating turbulence which is in accordance with the recent experiments. Also, a derivation for the 
two point third order structure function has been provided helping one to argue that even with slow 
rotation one gets, although dominated, a spectrum with the exponent -2.87, thereby hinting at the 
initiation of the two-dimensionalisation effect with rotation. Moreover, using the Gledzer-Ohkitani- 
Yamada (GOY) shell model, modifled for rotation, these signatures of two-dimensionalisation effect 
have been verified. 

PACS numbers: 47.27.1, 47.27. Jv, 47.32.Ef, 92.60.hk, 92.10.Lq 



I. INTRODUCTION 



Rotation, in the face of the discovery of two-dimensionalisation effect, has emerged as a parameter that can pro- 
gressively make a 3D turbulent flow look hke a quasi-2D or a 2D turbulent flow. The phrase 'look like' basically 
means that certain properties of 3D turbulence, such as wavenumber dependence of energy spectrum, direction of 
energy cascade etc., become such that they give impression that the flow is getting two-dimensionalised. In view of 
the fact that the dynamics of oceans, atmospheres, liquid planetary cores, fluid envelopes of stars and, other bodies of 
astrophysical and geophysical interest do require an understanding of inherent properties of turbulence in the rotating 
frame of reference, the problem of two-dimensionalisation is of central interest to any serious scientist; turbulence in 
rotating bodies is even of some industrial and engineering interest. 

In the steady non-turbulent flow, for low Rossby number [Ro = U/2Lil) and high Reynolds number {Re = UL/v), 
Taylor-Proudman theorem[l] argues that rotation two-dimensionalises the flow. This argument is often carelessly 
extended to turbulent flows to explain the rotation induced two-dimensionalisation arising therein. The two- 
dimensionalisation of the 3D turbulent flow in presence of rotation has begun to be understood as a subtle non-linear 
effect which is distinctly different from Taylor-Proudman effect. 

Cambon et al. [2] have showed that in the presence of rotation, the transfer of energy from small to high wavenumbers 
is inhibited; at the same time, the strong angular dependence of this effect leads to a draining of the spectral energy 
from the parallel to the normal wave vectors (w.r.t. the rotation axis) showing a trend towards two-dimensionalisation. 
Walefl^e[3] has used helical decomposition of the velocity field to study the nature of triad interactions in homogeneous 
turbulence and coupling it with the instability assumption predicted a transfer of energy toward wave vectors perpen- 
dicular to the rotation axis under rapid rotation. The helical decomposition turns out to be very handy to deal with 
rapidly rotating turbulent flow. In that case the linear eigensolutions of the problem, the so-called inertial waves, have 
the structure of helical modes. The assumption about the triadic transfers, coupled with resonance condition for non- 
linear interaction between inertial waves, show that there will be a tendency toward non-linear two-dimensionalisation 
of the flow. 
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Simulations by Smith et al.[4] speak volumes for the two-dimensionalisation effect. They show the coexistence of in- 
verse cascade (a typical feature of 2D turbulence) and forward cascade in forced rotating turbulence within a periodic 
box of small aspect ratio. In the simulations, the ratio of the mean rates of energy dissipated to the energy injected 
decreased almost linearly, for Ro less than a critical value, with decrease in Ro (increase in angular velocity |r2|). 
By the way, a very recent numerical study[5] shows similar transition from stratified to quasi-geostrophic turbulence, 
manifested by the emergence of an inverse cascade - a conclusion that agrees with that of Lindborg[6]. 
Although recent experiments by Baroud et al.[7, 8] and Morize et al.[9, 10] have shed some light on the two- 
dimensionalisation effect, the scaling of two-point statistics and energy spectrum in rotating turbulence remains 
a controversial topic. Zhou[ll] in analogy with MHD turbulence has proposed an energy spectrum E{k) ^ fc^^ for 
rapidly rotating 3D turbulent fluid (also see [12]) and this does seem to be validated by some experiments [7, 8] and nu- 
merical simulations[13-16]. But some experiments [9] do not tally with this proposed spectrum. They predict steeper 
than k~'^ spectrum and this again seem to be drawing some support from numerical results[17, 18] and analj^tical 
results found using wave turbulence theory [19, 20]. 

Unbiasedly speaking, if one wishes angular velocity to become a relevant parameter in constructing the energy spec- 
trum E{k), simple dimensional analysis would lead one to: 

Eik) oc O^s^fc-™ (1) 

where m is a real number, m should be restricted within the range 5/3 to 3 to keep the exponents of fl and s in 
relation (1) positive. The two limits m = 5/3 and m — 3 corresponds to isotropic homogeneous 3D turbulence and 2D 
turbulence respectively. The spectrum due to Zhou — E{k) ~ k~^ — is due an intermediate value of m = 2. So, as 
far as the present state of the literature on rotating turbulence goes, two-dimensionalisation of 3D turbulence would 
mean the dominance of a spectrum which goes towards E{k) ~ k~^ and which may choose to settle at E{k) ^ k~^ 
— an issue yet to be fully resolved. 

As a turbulent flow can be treated as the manifestation of a random velocity fleld, one actually hopes to unveil the 

statistical properties of the flow rather than every other detail of the flow. This means that we basically are after 
some probability distribution for the flow: Even today, this remains a tough nut to crack. However, the knowledge 
of the structure functions assists one to take the first step towards finding the distribution. The structure functions 
are experimentally measurable and hence are of extreme practical importance. The scaling relations of the structure 
functions, thus, are the lynchpins of turbulence theory although uncertainty lingers as to their general validity and 
the details of the derivation as far as the present status of research in turbulence is concerned. Therefore, naturally 
a lot of time and effort are spent by the scientists working in the field of turbulence to determine the exact forms for 
these functions and to study various phenomena in their light. 

In this article, we shall deal with the phenomenon of two-dimensionalisation of 3D incompressible high Reynold's 
number fiuid turbulence and try to see what can be said about it from the study of structure functions, especially S3 
(to be defined below). Basically, herein we shall comprehensively review the works [21-26] done in this direction by the 
author and observe how the method of calculating structure functions, developed by Kolmogorov and subsequently 
Landau, serves as the cornerstone for studying the two-dimensionalisation effect from the angle adopted by the author. 



II. 2D TURBULENCE 



It may be unanimously accepted that Kolmogorov's four-fifths law[27] is a landmark in the theory of turbulence 
because it is an exact non-trivial result. In three spatial dimensions, this law says that the two-point third order 
velocity correlation function behaves as: 



83 = 



|ti(x + 1) — v{x)^ .■ 



7el 



where e is the rate per unit mass at which energy is being transferred through the inertial range, v is the velocity field. 
The inertial range is the intermediate spatial region postulated by Kolmogorov where the large scale disturbances 
(flow maintaining mechanisms) and the molecular scale viscous dissipation play no part. This result is of such central 
significance that attempts are regularly made to understand it afresh and to extend it in other situations involving 
turbulence. There appears to be following two primary methods[28] of obtaining this result: 

1. The original Kolmogorov method put forward in details in the fluid dynamics text due to Landau and Lifshitz[29]. 
There is no external forcing in this approach, and the equality of dissipation rate and forcing rate for the energy 
is never enforced. 
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2. A ficld-thcorctic technique that invokes the so-called 'dissipation anomaly' in the high Reynold's number fluid 
turbulence. In this approach, there is an external forcing that maintains a steady state turbulence. 

The two approaches yield the same important results as they should. 

If one goes by the standard procedure given in the book by Frisch[30] to derive the form of the correlation function 
in d-D turbulence with the assumption of forward energy cascade, one would land up on[31]: 

where e is the mean rate of dissipation of energy per unit mass. This result is not quite true for the two-dimensional 
case since it gives for d = 2: ^3 = —{3/2)sl and not ^3 = (3/2)£Z. This is so because the calculation doesn't take into 
account the conservation of enstrophy in 2D turbulence which causes the reverse cascade of energy [32]. It might be 
noted that ^3 = {3/2)el for d = 2 is for the regime of scales larger than the forcing scale[21, 33]. 
Actually, if we consider the two-dimensional turbulence, then in the inviscid limit, we have two conserved quantities 
- energy and enstrophy. This gives rise to two fluxes with the enstrophy flnx occurring from the larger to the smaller 
spatial scales. The energy flux goes in the reverse direction. Recently, Bcrnard[33] and Lindborg[34] have used the 
above mentioned techniques to obtain the third order structure function for both the energy and the enstrophy cascade 
regions in forced 2D turbulence. We believe that the issue is important enough that a derivation of these results using 
the Kolmogorov-Landau approach should be useful. This is what we have attempted here and our results do come out 
in agreement with them. Besides, wc also have derived some other correlation functions which deal with vorticity fields 
in the inertial region and also some two-point second order correlation functions in the dissipative region following the 
arguments of Landau, thereby consolidating the equivalence between the two approaches mentioned in the beginning. 



A. Second order velocity correlation function 

It is a well-established fact that there exists a direct-cascade of enstrophy in 2D turbulence. One defines total 
enstrophy as F = ^ J^^ bP'Sp where w = d^Vy — dyVx is the vorticity in the Cartesian coordinates; v being the 

velocity field. As we shall consider incompressible fiuids only (V .v = 0), we shall take density to be unity and let p 
take over the task of representing position vector in 2D plane. The enstrophy flows through the inertial range and gets 
dissipated near dissipation scale. Using the antisymmetric symbol that has four components, viz. en = £22 = 
and £12 = —£21 = 1, one may deflne the mean rate of dissipation of enstrophy per unit mass as: 

T] = u{VijJ.Vu) (3) 

V = '^SraSe0{{drdjVa){d0dyV0)) (4) 

Here, angular brackets denote; an averaging proc;e(lure which averages over all possible positions of points 1 and 2 at a 
given instant of time and a given separation. Now, if vi and V2 represent the fluid velocities at the two neighbouring 
points at pi and p2 respectively, one may deflne rank two correlation tensor: 

Ba/3 = {{V2a " «la)(«2/3 " Vip)) (5) 

For simplicity, we shall take a rather idealised situation of turbulence flow which is homogeneous and isotropic on 
every scale a case achievable in practice in a vigorously-shaken-fluid left to itself. The component of the correlation 
tensor will obviously, then, be dependent on time, a fact which won't be shown explicitly in what follows. As the 
features of local turbulence is independent of averaged flow, the result derived below is applicable also to the local 
turbulence at a distance p much smaller than the fundamental scale. Isotropy and homogeneity suggests following 
general form for 

B^p = Ar{p)5c,0 + A2{p)py°0 (6) 

where Ai and A2 are functions of time and p. The Greek subscripts can take two values p and _L which respectively 
mean the component along the radial vector p and the component in the transverse direction. Einstein's summation 
convention will be used extensively. Also, 

P = P2-PI, Pa = PM, P°p = h p1 = (7) 
using which in the relation (6), one gets: 



= S±± (<5a/3 " + SppPaP^ (8) 
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One may break the relation (5) as 

Ba0 = {viaVl/}) + {V2aV20) - {viaV2p) - {V2aVl0) (9) 

and defining baff = (wiaW2/3), one may proceed, keeping in mind the isotropy and the homogeneity, to write 

Ba0 = {v^)Sal3 - 2ba0 (10) 

Again, having assumed incompressibihty, one may write: 

dpBa^p = (11) 
B^^ = pB'pp + Bpp (12) 

where the equation (8) has been used and prime (') is allowed to denote derivative w.r.t. p. Near the dissipation 
region the flow is regular and its velocity varies smoothly which allows to expand u in a scries of power of p. One 
must take v (? neglecting the higher powers (w ~ /9 is not taken because it leads to the contradictory result that 
?7 = as can be seen from the relation (4)). So, treating a as a proportionality constant, let Bpp = ap'^, which means 
B±± = 5ap^ (using equation (12)) and hence, 

1 5 

{viaV20) = -^{v'^)Saf3 - -^ap'^Sa/} + 2ap'^Pap/3 (13) 
{{dirdijVia){d20d2yV2i3)} ^ -72ad9rSa0 + 24:a600SaT + '24:a5aeS/3r (14) 

£Ta£ei3{{drd^Va){ded^vi3)) = -192a (15) 
- ^.. = -^ (16) 

In the equation (16), we have put p\~p2: for these relations are assumed to be valid for arbitrarily small p. While 
writing the relation (16), equation (4) has been recalled. This Bpp is the two-point second order correlation function 
for enstrophy cascade in dissipation range. 

B. Third order velocity correlation function 
Now, we shall focus thoroughly on the inertial range. Let's again define: 

Invoking homogeneity and isotropy once again along with the symmetry in the first pair of indices, one may write the 
most general form of the third rank Cartesian tensor for 6q;/3,7 as 

6a/3,7 = C(p)^„;3p; + D(p)(5^/3P° +^„^p^)+i^(p)p>^p; (17) 

where, C, D and F are functions of p. Yet again, incompressibihty dictates: 

^K,„ = = (18) 

2D^ 2D' 2D F' F ^ 

G ()„/3 H Oa/3 H Oap H irP<^Pli TP<^PP + -TPcPIS + -^PaP0 = (19) 

p p p^ P'^ p-' p-^ 

Putting a = /3 in equation (19) one gets: 

^ „ „ constant , , 

2C + 2D + F = = (20) 

where, it as been imposed that 6^/3,7 should remain finite for p = 0. Again, using equation (19), putting a ^ (3 and 
manipulating a bit one gets: 

D = -\{pC' + C) (21) 
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using which in relation (20), one arrives at the following expression for F: 

F = pC' -C (22) 

Defining 

= "^{^affa + &7/3,a + f'a7,/3) (23) 

and putting relations (21) and (22) in the equation (23) and using relation (17), one gets: 

B^p^ = ~2pC'{5^0P° + S^ppl + Sc^^pD + 6(pC' - C)pyip° (24) 
^ 53 = Bppp - -6C (25) 

which along with relations (21), (22) and (17) yields the following expression: 

6a/3,7 = - ^5o,0P" + ^ (P^S + ^3) {5^0Pl + 5o.-tPl) - \ (/^^3 ' ^s)^^^^" (^6) 

Navier-Stokes equation suggests: 

d 

d 

-Q-^'^W = -V2jd2^V2l3 - + vd2jd2jV2l3 (28) 

multiplying equations (27) and (28) with V213 and V\a respectively and adding subsequently, one gets the following 
after averaging the consequent result: 

d 

-dla{piV2l3) - d20{p2Vla)) 

+l^dijdij{ViaV2f3 + l^d2jd2j{viaV2fi) (29) 

Due to isotropy, the correlation function for the pressure and velocity, {{P1V2)), should have the form f{p)p/\p\. But 
since da{piV2a) = due to solenoidal velocity field, f{p)p/\p\ must have the form constant x (p/|p|^) that in turn 
must vanish to keep correlation functions finite even at p = 0. Thus, equation (29) can be written as: 

d 

Q^baf} = dj{ba'y,l3 + b(3j,a) + 2l'd^d^ba0 (30) 

For isotropic and homogeneous turbulence, the condition of incompressibility gives the easily derivable well-known 
result: 

Defining W = {CO1UI2) and noting that W = —dsdsbaa, we get from relations (30) and (31): 



dW _ 1 
'~dr ~ 2 



= -dsdsid^Bo^aj) - 2udsdsW (32) 



Again, if one defines 17 = {{lu2 — iu'i){LU2 — f^'i)) (which is not to be confused with the rotation rate discussed earlier), 
for homogeneous isotropic turbulence one may write il = 2(a;^) — 2W. So, equation (32) can be manipulated into the 
following: 

^ ^dsdsd^Baaj + vdsdsQ. - 2udsds{u'^) (33) 
dsdsd^B^aj = 47? (34) 

^ Baap = ^VP^ (35) 
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Here, we have assumed ^ to be relatively negligible and let ^ so that the terms proportional to v vanish. Also, 
we have recalled that = —J?- From equations (17) and (23), and the condition of incompressibility, it readily 

follows that B^^p = §-§^Bppp putting which in expression (35) and integrating subsequently (keeping in mind that 
Bppp shouldn't blow up at p = 0), we arrive at: 

Bppp = +lvp' (36) 

This is the one-eighth law for the unforced 2D incompressible turbulence proved using the Kolmogorov-Landau 
approach. 

Let us go back to equation (30). Using expressions (10) and (26), one can rewrite the equation as: 

As we are interested in the enstrophy cascade, the first term in the R.H.S. is zero due to homogeneity and the first 
term in the L.H.S. is zero because of energy remains conserved in 2D turbulence in the inviscid limit (and of course, 
it is the high Reynolds number regime that we are interested in); it cannot be dissipated at smaller scales. Also, as 
we are interested in the forward cascade which is dominated by enstrophy cascade, on the dimensional grounds in the 
inertial region Bpp (as it may depend only on r] and p) may be written as: 



dt 



-Bpp = Aj^p^ (38) 



where ^ is a numerical proportionality constant. Hence, using the relation (38), the equation (37) reduces to the 
following differential equation: 

1 a . 3 s V d f dBpp\ A 2 /on\ 
^SQ-^iP B,,,) = -g-p[P^)-2'^P (39) 

which when solved using expression (25) in the limit of infinite Reynolds number (i' — > 0), one gets 

Bppp = -^p' (40) 

Comparing it with the expression (36) for the two-point third order velocity correlation function for the isotropic and 
homogeneous 2D unforced turbulence in the inertial range of the forward enstrophy cascade, we determine the value 
of A to be -1/4. Therefore, relation (38) yields 

li^PP = -\vP^ (41) 

This, to the best of our knowledge, is yet another exact new result that has to be verified experimentally and 

numerically to test its validity. 

Suppose in the homogeneous isotropic fully-developed turbulence in two-dimensional space, energy is being supplied 
and the mean rate of injection of energy per unit mass is denoted by e. Let us concentrate on the inverse energy 
cascade. Then technically we have to proceed as before and on doing so one would re-arrive at the differential equation 
(37); only that now the arguments would differ. In the larger scales viscosity is not as significant and anyway we shall 
be interested in the infinite Reynolds number case which would mean that the last term in the R.H.S. of equation (37) 
would go to zero. One obviously would also set 5^(^^^) = £ and lets assume ■§fBpp w in the inverse cascade regime, 
justification of which can be sought from the fact that the ultimate result that is obtained has been experimentally 
and numerically verified. So, we are left with the following differential equation: 



1 d 

3 



^so-p(^'Bppp)=s (42) 



Bppp = +-ep (43) 

where in the last step the integration constant has been set to zero to prevent Bppp from blowing up at p = 0. This 
equation (43) is the expression for two-point third order velocity correlation function of the energy cascade in inertial 
the range. 
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C. Second order vorticity correlation function 

Recall that: 

W = {LO1LO2) (44) 
and, O = {{u>2 - 0Ji){ui2 - u>i)) (45) 

and that due to homogeneity, fl may be expressed as: 

n = 2{uj^} - 2W (46) 

In the dissipation range: t; ~ so, a; p and hence we may, choosing a proportionality constant b (say), assume: 

O = 6/9^ (47) 



Using relations (44), (46) and (47), one gets: 



But we know, 



So, relation (49) would yield: 



(6.1^.2) - (^') " (48) 

((5la'^)(52aW)) = 26 (49) 



T] = i^{{do,u;){daUj)) (50) 



rj = 2vb (51) 



where, we have put pi « p2 in the relation (49), for, being in the dissipation range, these relations are assumed to be 
valid for arbitrarily small p. Combining relations (47) and (51), one arrives at a experimentally verifiable result for 
two-point second order vorticity correlation function in the dissipation range: 

n = ±/ (52) 

D. Third order mixed correlation function 

Now, we wish to find two-point third order mixed correlation function in the inertial range of enstrophy cascade. 
We start by defining a two-point third order mixed correlation tensor in inertial range: 

% = {{V2(3 - Vii3){(jj2 - a;i)(w2 - wi)) (53) 
Q.p = 2Mf) + AWp (54) 

where, W/3 = (ui/3a'iW2) and Mg = {<jJ\uj\V2I}) ■ Due due to isotropy and homogeneity, we can write following form for 

= M{p)p"0 (55) 

d 

^ M/3 = {uJ\UJid2l3V2p) = (56) 

Op2t3 

- ^M(p) + ^=0 (57) 
, constant 

M(p) = = (58) 

P 

In the relation (56), we are assuming incompressibility and in writing the relation (58) we have taken into account 
the fact that should remain finite when p = 0. Relations (55) and (58) imply that: 

Mfs = (59) 
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using which in the relation (54), we get: 

= 4W/3 (60) 
Prom the equations (27) and (28), we may write respectively: 

d 

— wi = -t;i^ai^a;i + vdi^di^uii (61) 
d 

— U>2 = -V2jd2'yU>2 + vd2^d2'yU>2 (62) 

Multiplying equations (61) and (62) by 0J2 and ui respectively and adding subsequently, we get the following differential 
equation after averaging: 

^W = 2d0W0 + 2ud(}d(}W (63) 

where we have used the fact djs = —dijs = 8213 ■ Using relations (46) and (60) in the equation (63), one gets for the 
inertial range for the enstrophy cascade in homogeneous, isotropic and fully-developed freely decaying turbulence in 
two-dimensional space in the infinite Reynolds number limit {i.e., i/ — > 0) following differential equation: 

Q.p = -2r]p (65) 

In getting relation (65) from the equation (64), we have used the facts: = — rj and ~ as it may be 

supposed that the value of varies considerably with time only over an interval corresponding to the fundamental 
scale of turbulence and in relation to local turbulence the unperturbed flow may be regarded as steady which moan 
that for local turbulence one can afford to neglect ^fi in comparison with the enstrophy dissipation rate rj. This 
result (relation (65)) has gained importance by serving as the starting point in deriving various rigorous inequalities 
for short-distance scaling exponents in 2D incompressible turbulence [35]. 



III. QG TURBULENCE 

Quasi-geostrophic (QG) turbulence is a rather more realistic class of turbulent flow than the isotropic homogeneous 
3D turbulence. It can be seen in the large scale flows on oceans and atmosphere; thus having profound geophysical and 
astrophysical significance. QG turbulence [36] stands somewhere in between 2D and 3D turbulences. Thus, naturally it 
is very appealing candidate that deserves study if one is interested in the two-dimensionalisation effect. In the inviscid 
limit, besides total energy, QG flows enjoy the possession of yet another conserved quantity which is conserved at the 
horizontal projection of the particle motion. We shall call this pseudo-potential vorticity to distinguish it from the 
potential vorticity that is conserved at a particle in a homcntropic fluid. Defining pseudo-potential enstrophy as half 
the square of the pseudo-potential vorticity, one would say that like 2D turbulence there are two cascades — forward 
cascade of pseudo-potential vorticity and inverse cascade of energy — in QG turbulence which, however, is inherently 
three dimensional in nature. 

Recently, a paper [37] has calculated some structure functions in QG turbulence and has made illuminating revelation 
that isotropy in the sense of Charney[36] is useless in deriving the structure functions for QG turbulence. It has 
gone on to show that formulation of QG turbulence under the constraint of axisymmetry is productive. However, it 
criticized (though somewhat rightly) the ineffectiveness of use of tensorial quantities in the case of QG turbulence in 
deriving the results. Now, manipulating the tensorial quantities are at the heart of the derivation of many important 
two-point velocity correlation functions and other ones [29]. The technique is very intuitive and straightforward. It 
has, recently, also been thoroughly used to find out various correlation functions for 2D turbulence [21]. In this article, 
we shall closely (and trickily) follow the original Kolmogorov method put forward in details in the fluid dynamics 
text due to Landau and Lifshitz[29]; and repeated in the ref.-([21]), to derive structure functions in QG turbulence. 
The method has the extra advantage to being able to probe into the form for the two-point third order velocity 
correlation function in the forward pseudo-potential enstrophy cascade regime — this has remained uninvestigated 
earlier in ref.-([37]). 
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A. Third order mixed correlation function 



First of all, we shall briefly introduce the necessary equations (see ref.-([38]) for details). Let u be the three 
dimensional velocity field of the fluid in a frame rotating with constant angular velocity fl. The fluid body (such as 
ocean) is assumed to be of uniform density with free surface a.t z = ^{x,y,t). Suppose the bottom z = —H{x,y) is 
rigid. The shallow-water equations, then, are: 



dh 

— + W.{vh) = (66) 
and, ^+fxv = -gV^ (67) 

Here, h{x,y,t) = ^{x,y,t) + H{x,y), v = {ux,Uy), V = {dx,dy), -^^ = + v.V and v = v{x,y,t). f is Coriolis 
parameter that is Taylor-expanded to write / = /o + Py- Using the equations (66) and (67), one gets the relation: 



D 
Dt 



z.{curlu) + f 



h 



= (68) 



Let us assume: a) Rossby number Ro <C 1, b) Fractional changes in h are small, and c) /SL/fo <C 1 where L is the 
horizontal scale of the flow. Imposing these three assumptions on the shallow-water equations one can modify the 
relation (68) to yield 

^ + v.Vq = (69) 

where q = V^V + / ~ foi^/g^o + fo{Ho — H)/Hq {tp being g^/fo) may be called pseudo-potential vorticity. Under 
the same assumptions, for QG flow, one also has the condition: 

V.v = (70) 

Now, the trick is to select an arbitrary two-dimensional plane in the QG turbulent flow such that the plane's normal is 
along / and impose the properties of homogeneity and isotropy in the plane only. By the way, one must keep in mind 
that the so-called fundamental scale of 3D turbulence has its analogy as the horizontal length scale L for the case of 
QG turbulence. The correlation functions to be derived for the forward cascade in this article are valid in the range 
(which we shall call inertial range) that is much smaller than L but quite larger than the scale at which the dissipation 
is eff'ective. Whereas, the structure function to be derived for the inverse cascade is valid in the range whose scale is 
larger than the scale at which energy had been fed in. As we shall consider fluid bodies of uniform density only, we 
shall take density to be unity and let p, as usual in this article, take over the task of representing position vector in 
the 2D plane. The Greek subscripts used herein can take two values p and _L which respectively mean the component 
along the radial vector p and the component in the transverse direction. Whenever we shall use the Latin subscript 
{e.g., a), it should mean that it can take one more value apart from the ones mentioned above: the third value 'z' 
would signify the vertical direction. As before, Einstein's summation convention will be used extensively. Also, 

P = P2-Pl, p°^ = pa/\p\, P°p = l, P1=0 (71) 

Now, if vi and V2 represent the horizontal fluid velocities at the two neighbouring points at pi and p2 respectively 
then with similar meaning for qi and q2, one may define, just for the sake of notational convenience: 

K = {qiq2) (72) 
and, Q = {{q2 - qi){q2 - qi)) (73) 

The angular brackets denote an averaging procedure which averages over all possible positions of points 1 and 2 at a 
given instant of time and a given separation. Due to homogeneity, Q may be re-expressed as: 

Q = 2{q^) - 2K (74) 

For simplicity, we shall take a rather idealised situation of QG turbulence which is homogeneous and isotropic on 

every scale in the plane. For the unforced case, the component of the correlation tensor will obviously be dependent 
on time, a fact which we won't be showing explicitly in what follows. As the features of local QG turbulence should 
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be independent of averaged flow, the result derived below is applicable also to the local turbulence in the plane at 
scale p much smaller than the fundamental scale. 

Again, we define a two-point third order mixed correlation tensor in inertial range: 

Ql3 = {{V20 - Vip){q2 - qi){q2 - qi)) (75) 
^ Q0 = 4K0 + 2L0 (76) 

where just to reduce the effort of writing, we have defined: 

Ki3 = {viuqiq^) (77) 
and, Lj3 = (gigi^/s) (78) 

Obviously, isotropy, homogeneity and the condition (70) compels L/3 to vanish. Hence, equation (76) reduces to: 

(3/3 = (79) 

From the equation (69), we may write for the points 1 and 2 respectively: 

d 

—qi = -vi^di^qi (80) 
d 

—q2 = -V2jd2jq2 (81) 

Multiplying equations (80) and (81) by q2 and qi respectively and averaging subsequently after adding, we get the 
following differential equation: 

^=2d0K0 (82) 

Using relations (74) and (79) in the equation (82), one gets for the inertial range for the pseudo-potential enstrophy 
cascade in homogeneous and isotropic QG turbulence (forced at an intermediate scale or unforced) in inviscid limit 
the following differential equation: 

=> Qp = -2egp (84) 
In getting relation (84) from the equation (83), we have assumed the following: 

1. Qp does not blow up at p = 0. This sets the integration constant as zero. 

2. ^-gfiq^) = —£q, i-e., there exists a pseudo-potential enstrophy sink at small scales due to some dissipative force 
such as viscosity and £q is the finite and constant dissipation rate of the mean pseudo-potential enstrophy. 

3. ~ due to quasi-stationarity. It may be supposed that the value of Q varies considerably with time only 
over an interval corresponding to the fundamental scale of turbulence and in relation to local turbulence the 
unperturbed flow may be regarded as steady which mean that for local turbulence one can afford to neglect 

in comparison with the pseudo-potential enstrophy dissipation rate Sq. 

B. Third order velocity correlation function 

Having explored the form for two-point third order mixed correlation function in the preceding discussion, we now 
proceed to find the scaling for the two-point third order velocity correlation function. For this motive, one may define 
a rank two correlation tensor: 

Ba0 = {{V2a - Vla){V20 - Vip)) (85) 

Isotropy and homogeneity in the plane suggests following general form for 

Ba^fS ^ Aiip)6af3 + A2{p)p°„p°0 (86) 
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where Ai and A2 are functions of time and p. Making use of the relations (71) in the equation (86), one gets: 

Ba0 = Bi±(<5a/3 - P>^) + -Bp,P>^ (87) 

One may expand the R.H.S. of the relation (85) and defining ba^ = {viaV2i3), one may proceed, keeping in mind the 
isotropy and the homogeneity, to arrive at: 

Ba/3 = {v'')Sa0 - 2baf3 (88) 

Let's concentrate on the following statistically averaged quantity that will prove to be of crucial importance for 
deriving the desired results: 

ba/3,-f = {vioiVll3V2-y) 

Invoking homogeneity and isotropy in the plane once again along with the symmetry in the first pair of indices, one 
may write the most general form of the third rank Cartesian tensor for this case as 

= C(p)^„;3p; + £>(p)(^^^p° +^„^p^)+i^(p)p>^p; (89) 

where, C, D and F are functions of p. Imposing the condition (70) on the expression (89), one can get (in the same 
way as done earlier for the 2D case) the following relations: 

D = -\{pC' + C) (90) 
and, F = pC -C (91) 

Here, prime (') denotes derivative w.r.t. p. Defining 

= '^{ba(3,j + ^7/3,a + ^'a7,/3) (92) 

and putting relations (90) and (91) in the equation (92) and using relation (89), one gets: 

= -2pC"(^„;3p^ + 5^0pI + 5^^pI) + 6(pC' - C)plplp° (93) 
which along with relations (89), (90) and (91) yields the following expression: 

6a/3,7 = + + B,,,){5^ppi + 5^,pi) - i(pi?;,^ - i?,,,)p>^p; (94) 

The equation (67) suggests: 

d 
d 

01""^ = -V2-fd2jV2l3 + /ioeo/37^27 - gd20^2 (96) 
multiplying equations (95) and (96) with V20 and via respectively and adding subsequently, one gets the following: 

-g^{viaV20) = -dl^{vijViaV2l3) - d2^{V2^ViaV2l3) 

+ {flaVl^V2f3) + eal3j{f2aV2'yVla) 

-gdla{^lV2f3) - gd2f3{^2Vla) (97) 

Due to isotropy, the correlation function {£,iV2) should have the form f{p)p/\p\. This / should not be confused with 
the Coriolis parameter. But since, 5a(^ii;2a) = owing to the relation (70), f{p)p/\p\ must have the form kp/\p\'^, 
where fc is a constant. Now, k must vanish to keep correlation functions finite even at p = 0. Thus, equation (97) can 
be written as: 

d 

— 6a/3 = d.y{baf,l3 + bi3f,a) + fo^zaibjff + fo^zfj-ybai (98) 
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Here we have used the approximation: / = foz. Using equations (88) and (94), one can rewrite equation (98) as: 

Note that the terms containing the Levi-Civita symbol vanish by virtue of the joint effect of the expressions (87) and 
(88), and the antisymmetry property of Levi-Civita symbol. As we are interested in the pseudo-potential enstrophy 
cascade, the first term in the L.H.S. is zero bccaiisc of energy remains conserved in QG turbulence in tlio inviscid hmit: 
it cannot be dissipated at smaller scales. Also, as we are interested in the forward cascade which is dominated by 
pseudo-potential enstrophy cascade, on the dimensional grounds in the inertial range Bpp (if it is assumed to depend 
only on Sg and p) may be written as: 

^Bpp = Tsy (100) 

where F is a numerical proportionality constant. Hence, using the relation (100), the equation (99) reduces to the 
following differential equation: 

= (101) 
which when solved, imposing finiteness of Bppp for p = 0, gives: 

Bppp = -^p' (102) 

The relation (102) is the expression for the two-point third order correlation function in the isotropic and homogeneous 
plane of QG turbulence (forced or unforced) in the range of the forward cascade where there is no overlapping with 
energy cascade. Since F has not been determined, one must confess that the equation (102) is just a scaling law at 
this stage. 

Now suppose the fluid body is being forced at small scales i.e., energy is being supplied and the mean rate of injection 
of energy per unit mass is denoted by £„ (assumed finite and constant). Let us focus on the inverse energy cascade. 
Then technically we have to proceed just as before to finally arrive at the differential equation (99). One obviously 
would set \^{v'^) = invoking the hypothesis [36] that there should be equipartition of energy between potential 
energy and the energy content in each of the two horizontal velocity components in the plane. This equipartition 
had been proposed in view of the assumption that at sufficiently small scales the interaction of the mean flow with 
the eddies (and thus the eddy-cnergics) diminishes; as a result, for increasingly smaller vertical and horizontal scales 
the energies will tend to become homogeneous and equally distributed among the perturbations. By the way, the 
concept of equipartition of energy is very old and wide-spread in the literature of statistical mechanics. Historically, 
equilibrium statistical mechanics had been used to justify many aspects of turbulence e.g., the dual cascades in 2D 
turbulence etc. A detailed discussion may be foimd in the books by Chorin[39] and Lim et a/. [40]. Now, lets also 
assume that §iBpp « in the inverse cascade regime supposing the forced QG turbulence to be in the state of 
quasi-stationarity. So we are left with the following differential equation: 

^ Bppp = +SuP (104) 

where in the last step the integration constant has been set to zero to prevent Bppp from blowing up at p = 0. 
The expression (104) is the expression for the two-point third order velocity correlation function in the isotropic and 
homogeneous plane of forced QG turbulence for the inverse energy cascade. 

The fact that the structure functions for the inherently three-dimensional QG turbulence are more like that for the 
2D turbulence than that for the 3D turbulence speaks volumes for the importance of study of third order structure 

functions for demystifying the two-dimensionalisation effect of the 3D turbulent fluid due to rapid rotation. This 
serves as the motivation for jumping into the subject of rotating flows and to attempt finding the form of 6*3 therein. 



IV. ROTATING TURBULENCE 



All the studies on the two-dimensionalisation effect are mainly for low Ro high Re limit while the high Ro and 
high Re limit has been rather less ventured in relation to the two-dimensionalisation effect of turbulence, although 
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the second case, we believe, should be analytically more tractable. If, using calculations of structure functions, in the 
limit of high Ro and high i?e, one wishes to see whether a trend towards two-dimensionalisation of 3D homogeneous 
isotropic turbulence occurs or not, then basically one would have to check (a) if 6*3 — — (4/5)e/ at small scales for 
3D turbulence shows a tilt towards = {3/2)el at large scales for the 2D turbulence and (b) if the forward energy 
cascade is depleted at the smaller scales. As we shall show, in the lowest order calculation this is what one may get, 
hinting at the initiation of the effect of two-dimensionalisation of 3D turbulence owing to the small anisotropy induced 
by slow rotation. 

A. Relevant scales 

Let us look in to the various length scales that have to be taken into consideration while talking about a homogeneous 
rotating turbulence which basically satisfies following version of Navier-Stoke's equation: 

^ + (w.v) V = -^VP -dx (nxxj -2nxv + z/V^u + / (105) 

In this context / is external force and fl is angular velocity. Various parameters to be considered are: v (kinematic 
viscosity), e (finite mean rate of dissipation of energy per unit mass), O (angular velocity) and Iq (integral scale 
which typically is the system-size). The three important time-scales involved in the system arc: ti - e-i/'^/^/a (eddy- 
turnover time or circulation time for the eddy of scale /; I < Iq), ta ~ and td ~ l^/i^ (diffusion time scale). It is 
well-known that a length scale Iq = ^/{e/iW) is what responsible for the estimation of the anisotropy introduced by 
the rotation. The competition between the time-scales t; and gives rise to what is known as dissipation length scale 
Id, defined as Id = (i^'^/e)^/** and a similar competition between the time-scales td and tn allows us to define a length 
scale lad = \/ i'^/^)- Now, lets look at the typical scenario when Ro is moderate. The four vital length scales are 
typically arranged according to the order : Iq > Iq > l^d > Id- Thus, the regime Iq > I > Iq is the regime where effect 
of rotation is important and anisotropy reigns. The scales I G {la, Id) may be considered to have isotropy, though 
to be precise, probably Id here should be replaced by Ind since rotation seems to be bringing the effect of viscosity 
to rather larger length scales. So, now what happens when the Ro is decreased by increasing the angular velocity is 
interesting. Both the scales Iq and kid rush towards the dissipation length scale, thereby increasing the anisotropic 
regime and at the angular velocity fl — Qa = \/ i^/^) one has Iq = lad = Id and the turbulence is fully anisotropic. 
Strictly speaking, even a small rotation introduces anisotropy (however small) at all scales and the isotropic regime 
does have a degree of anisotropy in it as we shall see shortly. In the fully anisotropic limit, i.e. for f2 = f2o, one expects 
full decoupling of the plane perpendicular to the rotation axis from the direction of the rotation axis. However, even 
in the partially anisotropic limit {e.g. when we have slow rotation imparted on the turbulent fluid), Iz should still be 
given a special status for being in the direction of the rotation axis, by which we mean that the structure functions 
should no longer depend on I but rather on Iz and l±_ (where P = l'^ + and jilj = flz). 

We shall see how this decoupling sets in, in the limit of low angular velocity and try to study in that very limit, the 
two-point third order structure function in the first approximation and see how the effect of two-dimensionalisation 
is all set to sneak in with the switching on of rotation. 

B. S3 for small Q 

Let us start with low 0-limit. With this statement we mean, as discussed in the previous subsection, 17 <C rio- 

So, the entire fiuid may still be treated as isotropic but as rotation should play a role, we assume that (wiWjW^) 
(where angular brackets mean ensemble average and Vi = Vi {x, t) is the i-th component of velocity and similarly, 
vl = Vi{x + l,t)) should depend on Q as well. Q would take care of the mild anisotropy. Since, physically speaking, 
S3, should not depend on which way the rotation axis is and since we are interested in low values of 17, we shall let 
{viVjv'j^ depend only on the terms quadratic in 12 and not bother about higher order terms in 12. As a result, we 
write the following most general tensorial form for {viVjv'^,): 

bij,k = {viVjv'k) (106) 
= C{l)SijF, + D{l){5ikl] + SjkO + F{l)imil 

+G(0[(e^™fe^° + ^3mki°)i°J + H{i)n,n,ii 

+I{l)[{eimk^j + ejmki^i)nm] + K{l){niflkl° + 17,17fc/- ) (107) 
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where 1° is the i-th component of the unit vector along /. We have assumed that the coefficients arc dependent only 
on / and it is the ^ which is taking care of the mild anisotropy which the turbulent fluid might have. We must 
accept that the assumption of letting coefficients depend only on / is rather crude in the light of the complex forms 
that the two-point tensors in a fully anisotropic turbulence flow take[41]. The justification, and hence solace, for 
the assumption can be drawn from the fact that very simple revealing results matching with recent experiments are 
arrived at in the long run. As we are considering incompressible fluid, we must have: 

d'khjM = (108) 

which when applied to relation (107), yields relationships between various coefficients. As usual Einstein summation 
convention has been extensively followed in these calculations unless otherwise specified. Using relations (107) and 
(108), one lands up in the end on the following: 

Bijk = {{v'i-Vi){v'^-Vj){v'^-Vk)) (109) 

= 2(6y,fe + hjk,i + hi,]) (110) 
= ^2{IC' + C){5,,ll + 5,kl° + 5jkl°) + Q{IC' - C)l°l°ll 

+'iJl{Vl,Vljll + VL,Vtkl° + n^Vlkl") (111) 

Here, in expression (111), prime (""') denotes derivative w.r.t. / and J is a constant which, curiously enough, is of the 
same dimension [L^T~^] as that of the kinematic viscosity. Now we can see that using the relation (111), two-point 
third order structure function [S^) can be extracted from Bijk in the following way: 



Ss{l) ^ {{5v\\{l)f) 



^v{x + 1) - ■?(x)| .J 



(112) 



S^{1) = ((K - Vi)mv'^ - v^)^][{v'^ - v,,)ll]) (113) 
53(0 = Biikl°lTk (114) 
Ss{l) = -12C+-^{n.lf (115) 



where we have used relation (111). One may define physical space energy flux {s{l)) as: 

6(1) ^ -^Wi.{\Sv{l)\'Sv{l)) (116) 

^ e{l)=lC" + 7C' + ^ + 3Jn^ + ^{nIf (117) 

To get relation (117), we have again made use of the relation (111). The energy flux through the wave number K 
(Ilif) for the isotropic homogeneous turbulence may be calculated to be: 

IIk = - r dl'^^^il + ldi)s{l) (118) 
7^ Jo « 

Now if one makes the standard assumption (often made during the derivation of S3) that as Re 00, the mean 
energy dissipation per unit mass £{i') tends to a positive finite value {i.e., lim,^^o £(^') = £ > 0), then lim,^_>o = £ 
in the inertial regime. Therefore, in the inertial range, putting x = Kl, one has 

„ 2 r°° , sinix) , f x\ 

n. = -y^ do^^f{-^)=e (119) 

where, 

f[^)=f{l) = {l + ldi)e{l) (120) 
For small I (large K), the integral in relation (119) yields 



f{l)^e 



(121) 
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Now using relations (117), (120) and (121), wc form a differential equation which when solved, keeping in mind that 
S3 should not blow up at ? = 0, one gets following form for 5*3 in slowly rotating homogeneous turbulent fluid. 

Ssil) = -U + ^Jlp" + Ti^kllf] (122) 
o 

One may note from the relation (122) that how magically O has brought up the anisotropic effects even for small O 
though for the entire calculation we followed the procedure meant for the homogeneous isotropic turbulence. Thus, 
the form for S3 is plausible. 

One may ask: Does the effect of two-dimensionalisation shows up in the relation (122)? As one may note from the 

relation (122) this is quite a possibility but the only catch being that J should be positive an issue which wc have 
not been able to resolve. If J is positive, it means if we increase Q. the value of S3 would distort away from the usual 
— (4/5)£Z for the non-rotating case to more positive values. This apparently shows that the effective value of e is 
decreased depicting that the forward energy transfer is depleted which is in keeping with what is exp(X't{^d and hence 
the tendency of the rotating 3D turbulence to show the effect the two-dimensionalisation is being liigliliglited. That 
the sign of J should be positive is a question remains to be addressed. 

By the way, the relation (122) also suggests that the coefficients in the tensorial form for hij^k should have dependence 
on Iz and l±_ separately effecting a mild decoupling of directions. So taking hint from it, we proceed to rewrite hij^k 
for slowly rotating 3D turbulent fluid but now introducing anisotropy directly into the coefficients and not letting Q. 
take care of anisotropy explicitly. Of course, the coefficients will now depend on Q,. 

For completely isotropic homogeneous turbulence, one would write following general form (relation (123)) for hij^k 
which is made up of Kronecker delta and components of the unit vectors 

bij,k = C{l)5ijll + D{l){5ikl°j + 5jkl°) + F{l)ltl°j^ (123) 

The expression is symmetric in i and j and the coefficients are dependent on I only. As discussed earlier, with rotation 
coming into effect, anisotropy comes into effect. If this effects in the possible decoupling (even if partial) of the 
direction along the rotation axis (which we shall take along the z-axis), then mathematically we may introduce this 
effect by modifying the form (123) of bij,k to the following: 

hij,k = C{1, n)6,jll + D{1, h,n)id,kl° + Sjkli) + F{1, h, (124) 

If one uses the incompressibility condition (relation (108)), one gets: 

D=i{-C'-^)-C (125) 
and D = (126) 

where dot represents the derivative w.r.t. Iz and prime, as before, the derivative w.r.t. I. Using equation (125) in the 
equation (126), one land up on: 

Clz + IC' + 3C = (127) 

=> C = Y,Anl~"~^lz (128) 

n 

^ C ^ for n G (-00, -2] n [0, 00) (129) 
C = D = F = (130) 

In arriving at the result (130), we have taken care of the fact that C can not be allowed to blow up for either for 
^2 = or for I = 0. Thus, relation (124) vanishes trivially. So, we are left with the following choice: 

b^,,k = C{hJz,^)Sijll + D{hJz,mS^kl] + Sjkl-) + F{l^,lz,n)l°l°ll (131) 

Using equations (108) and (131), we arrive at following relationship between the coefficients: 

D = J-^C--C-C (132) 
2 2 ^ ^ 
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Here tilde and dot define derivatives w.r.t. l± and Iz respectively. Proceeding monotonously as before we get 

Bijk = '^{bij.k + bjk,i + bkij) (134) 
= -2{l^C + hC + C){5ijll + 6ikl] + Sjklt) + QFliljlk (135) 

And hence, 

^3 = BijklTjik = - + hC + C)] (136) 
The definition for the physical space energy flux {s{l)) has to be obviously modified. Natural choice would be: 



{\5v{r)\^Sv{l)) = BuJ"J-j^ + BuJ°J-^ 



(137) 



where a takes two values: x and y only. Now, using relations (116), (133), (135) and (137) and performing tedious 
algebra one gets: 



I 



{311 + ^iJz + ^liQc 



+ {3lil, + 6llll + 31^11 + Ql\ll + + 3«f) (5 

+ {2,1111 + (SI Jl + ?,iiHl)d 
+ (5/i + + 23i\il + mill + mjl + 6il) b 
+ {-ll\l^ + - + + 6Zf + i8z;^;^) c 
+ {-\2i\ - siliz - + mjl + mt + si^^il) c 

+ {-m\l, - iSllll - 39l±ll - 171*) c 
+ {-All - 8?i«. - I2l±ll) C] 



(138) 



The energy flux (11;^) through the wave number K for the homogeneous (not necessarily isotropic) turbulence may 
be shown to be: 



^3,sin{Kl) 



/M3 ^ 

Using cylindrical polar coordinates we reduce the relation (139) to: 



-(4 



1 



l_\_dl^dlz 



ii\n{Kl) 
I 



Ix d h d 1 



e{l) 



(139) 



(140) 



Now, we introduce the variables y = Kl± and z = Kl^ in relation (140) to get: 

sin(y^ + z^)^ 



IIk = - / 

Jz=-oo Jy=0 



dydz- 



y"^ + 



f(]L 1.] 

\K' kJ 



Now, let's probe small I behaviour. Because J^_^ J^q dydz[sm{y'^ + z^Y^"^]/ (y^ + z"^) = we have 



/(/±,y«- 



(141) 



(142) 



Obviously, e has the meaning of finite positive mean rate of dissipation of energy per unit mass. Using the expressions 
(138) and (142), we look for the Iz = hmit. One then has the result: 



ol± 



(3ZiC + U^C - 12(7 - 4^) = - — 

tX TT 



8£, 



3l\C + 14Z1C - 211C - 16Z_lC = Z_l 

TT 



c 



el± 



7-V97 7+^97 . ^, . 

+ A2lZ^ + A^l^ = + AJ^ ' ) + ^ 



(143) 

(144) 
(145) 
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Relations (133), (136) and (145) together yield following expression for S3: 



.^3^=0 = — Sl^+Ai 
TT 



7 + V97 \ / 1 + V97 



12 



7+V97 



Ssli^^o = el±+Al 



7+V97 

6 

_L 



(146) 



(147) 



where, ^ is a constant which for obvious reason depends on and s. Using dimensional arguments and introducing 
a non-dimensional constant c, we may set 



A = Cil * £ 12 



From relations (147) and (148), we may write finally 



^ , 6 , , ^ l + V^ 11-V97 I±vgZ 

TT 



(148) 



(149) 



This (relation (149)) is the final form for two-point third order structure function in the plane whose normal is parallel 
to the rotation axis for slowly rotating homogeneous 3D turbulence. 



C. Energy spectrum for small Q 

If we for the time being forget about the issue of anomalous scaling, then a mere inspection of the relation (149) 
from the point of view of dimensional analysis would tell that in the directions perpendicular to the axis of rotation, 
there are two possible energy spectrums viz. 

E{k) ~ fc-t (150) 
and, E{k) ~ k'^^^ (151) 

which are respectively due to the first term and the second term in the R.H.S. of the relation (149). It is very 

interesting to note that the exponent of k in the relation (151), i.e. — (16 + \/97)/9, equals —2.87 which is in between 
—3 (for 2D turbulence) and —2 (for rapidly rotating 3D turbulence as proposed by Zhou). Obviously, the spectrum 
(150) will be dominant compared to the spectrum (151). But as the CI is increased (of course, remaining within a 
range so that the anisotropy is not strong enough to breakdown the arguments used to calculate the S3 of the relation 
(149)), the spectrum (151) becomes more and more prominent; thereby two-dimensionalisation of the 3D homogeneous 
turbulent fluid is initiated which then carries over to high rotation regime as is being extensively studied now-a-days. 
This signature of two-dimensionalisation is, of course, in agreement with what present literature on turbulence hails 
as the two-dimensionalisation of turbulence. Thus, the third order structure function has proved to be very handy in 
studying this effect because the phenomenon of two-dimensionalisation is reflected as a change in the scaling law of the 
third order structure function. We pause here for a moment and ponder upon the signatures of two-dimensionalisation 
effect in a rather more intuitive, though a bit non-rigorous, way. 



D. The signatures: an intuitive picture 



Let us first concentrate on why at all there should be an inverse cascade of energy. Inverse cascade of energy is a 
trademark of 2D turbulence where a second conserved quantity — enstrophy — besides energy plays the defining role 
behind it. One might be tempted to search for this conserved quantity in the case of rapidly rotating 3D turbulence, 
for, there in the limit of infinite rotation the axes of all the vortices are expected to point up towards the direction 
of angular velocity. Hence, looking at the every section perpendicular to the axis one might tend to think that 2D 
turbulence is being shown by each transverse section. This obviously is not a correct inference because of the non-zero 
axial velocity may depend on the coordinates on the plane. Searching for the enstrophy conservation seems to be a 
dead end as far as explaining the inverse cascade in rapidly rotating turbulence is concerned. In such an unfortunate 
scenario, helicity (defined as / v.Lod^r) which remains conserved in a 3D inviscid unforced fiow comes to our rescue. 
It has been long known that helicity is introduced into a rotating turbulent flow [42]. Kraichnan[43] argued that both 
the helicity and energy cascade in 3D turbulence would proceed from lower to higher wave numbers and went on to 
remark that forward helicity cascade would pose a hindrance for the energy cascade — a fact validated by numerical 
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FIG. 1: A knotted vortex tube. When it is stretched the tube thins out to create smaller eddies but the entire structure 
occupies a larger volume. 

simulations [44, 45]. He also showed that in presence of helicity two-way cascade is possible. Lets see topologically why 
this should be so. It is well known that a knotted vortex tube is capable of introducing helicity in fluid[46]. Consider 
a knotted vortex tube (see Fig-1) in a turbulent flow. Due to the vortex stretching phenomenon in turbulence, the 
vortex line stretches and as a result owing to the assumed incompressibility of the fluid the tube thins out keeping 
the volume inside it preserved and smaller scales are created; in a sense, this is what is meant by the flow of energy to 
the smaller scales. But now this also means that the "scale" of the knotted structure would in general increase i.e., 
the knot would now reach out to farther regions in the fluid. Evidently, if we wanted this scale to reduce, we must let 
the stretched knotted tube fold in such a way so that the scale becomes smaller; such a neat arrangement seems to be 
a far cry in a turbulent flow which is inherently chaotically random causing the separation of two nearby particles of 
fluid on an average. Thus, as the degree of knottedness measures helicity, the aforementioned argument suggests that 
if one forces energy to go to smaller scales, helicity would tend to go to larger scale and vice-versa. This topological 
argument gives an intuitive way of comprehending how the forward helicity cascade can inhibit the forward cascade 
of energy. The point is that in presence of forward helicity cascade, reverse cascade of energy is not impossible. 
Waleffe[47], with the help of detailed helicity conservation by each triad, showed that helicity indeed affects the 
turbulence dynamics even in isotropic turbulence; this is a kind of catalytic effect. One can thus take inspiration to 
make the argument in the previous paragraph more concrete by playing around with a simplified triad using logic 
in the line suggested by Fjortoft's theorem[48] in 2D turbulence. Let the helicity spectrum be H{k) and the energy 
spectrum be E{k). It may be shown that 

\H{k)\ < kE{k) (152) 

Consider 3D Euler equation in Fourier space truncated in order to retain only three parallel wave vectors fci, 

and fca and suppose it is possible for these three particular wave vectors to be such that \H{k)\ = nkE{k), where 
n is a positive number lesser than 1 to be in consistence with the relation (152). Assume = 2A;i and ks = 3ki. 
Conservation of energy and helicity imply that between two times ti and t2, the variation SEi = E{ki, 12) — E{ki,t\) 
satisfies two constraints 

6E1 + SE2 +SE3 = (153) 
and, nkiSEi + nk25E2 + nk^SEs = (154) 

solving which in terms of SE2, we get: 

SE2 

SEi = 5E3 = ^ (155) 

and, nkiSEi = —^k2SE2; nk^SE^ = — ^fc2<5£^2 (156) 

If one assumes that the wave vector k2 is losing energy, i.e. 6E2 < 0, then the results (155) and (156) show that as 

more helicity goes into the higher wavcnumbcr. the energy is equally transferred to both the lower and the higher 
wave numbers suggesting a possibility of the coexistence of reverse and forward energy cascades. 
Now let us come to the point. In the case of 3D isotropic and homogeneous turbulence rotation can input helicity in 
it when there is a mean flow in the inertial frame and this value of input helicity increases with the increase in angular 
velocity. Experiments on rotating turbulence invariably introduce helicity. As the angular velocity is increased the 
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helicity increases enough to inhibit the energy cascade appreciably so that a reverse cascade is seen. This consistently 
explains the reason behind the existence of the reverse energy cascade in a rapidly rotating turbulent flow. Hence, 
the argued existence of a direct helicity cascade in such experiments turns out to be an interesting (however not 
rigourously proven) assumption. 

As discussed earlier, the next important signature of the two-dimensionalisation of turbulence that remains to be 
pondered upon is the exponent of the wave vector in the energy spectrum relation. To be precise, if one wishes 
angular velocity becomes a relevant parameter in the energy spectrum E{k), simple dimensional analysis would give: 

E{k) oc n^e^fc-" (157) 

where m is a real number, m must be restricted within the range 5/3 to 3 to keep the exponents of fl and s in 

relation (157) non-negative. The two limits m = 5/3 and m = 3 corresponds to isotropic homogeneous 3D turbulence 
and 2D turbulence respectively. The spectrum due to Zhou - E{k) ~ — is due to an intermediate value of 
m = 2. So, as far as the present state of the literature on rotating turbulence is concerned, two-dimensionalisation of 
3D turbulence would mean the dominance of a spectrum which goes towards E{k) ~ k~^ and which may choose to 
settle at E{k) ~ k'"^. 

Lets give a twist to the tale. In general, the energy spectrum[42] in the inertial range will be determined by both 
the helicity cascade and the energy cascade which simply means that the energy spectrum from the dimensional 
arguments should be written as 

E{k) oc e3-'"/i'"-tfc-'" (158) 

where h is the rate of helicity dissipation per unit mass. Demanding positivity of the exponents of e and h, one fixes 
the possible values for k within the closed range [5/3,7/3], imposing which on the arguments given in the previous 
paragraph, one can easily propound the range 

7 

2 < m < - (159) 

for the rapidly rotating 3D turbulent flow. Direct experiments [9] by Morize et al. have found energy spectrum for 
rapidly rotating turbulence going as k^"^'^ which is as predicted by the relation (159). 

One may note that the scaling exponent derived as expression (151) has not fallen into the more strict range 
[—7/3,-5/3] obviously because 17 is too low and may be because to maintain isotropy to a certain extent for the 
sake of hiccup-free calculations we have chosen not to include terms involving e^j]^ in the relation (131) which could 
grab the effect of helicity explicitly; thereby again showcasing the need for the helicity to be effective to give the right 
exponent for the rotating turbulence. 



E. Yet another signature 



Having explained the two signatures of the two-dimensionalisation effect, we search for another possible signature of 
the effect. The advection of a passive scalar 9 may serve the purpose since the Yaglom's law[49] in d-D incompressible 
turbulent fluid may be written as {5v\\{59)'^) = —{4:/d)eel, where ee = K{{^l^9) {di.9)) = —dt{9^) and k being the 
diffusivity. This law distinguishes between a 2D and a 3D turbulence and hence it is worth getting a form for it 
for a rotating 3D turbulence and find if in a plane perpendicular to the rotation axis it reduces to the form for 2D 
turbulence and thereby bringing in the effect of two-dimensionalisation. Since we have witnessed earlier that small f2 
could bring in anisotropy in the otherwise isotropic scales, one would look out for the effect of small on the passive 
scalar which follows the equation: 

fi9 f) 

— + y.{v9) = kS/H - eijk^j — {xk9) (160) 

If one goes by the procedure given in the reference [50] to find out a value for (i5f||(i5^)^) for small / in this case 
assuming very small fl (and hence isotropy), one arrives back at the Yaglom's law. We can however land up on a 
very neat experimentally and numerically verifiable correlation which can serve the purpose of a signature of two- 
dimensionalisation if we treat equation (160) anisotropically as follows. 

Defining I = x' — x and di^ = Vj = 9- = —di, one can manipulate the equation (160) to get: 



dt{{69f) + Vi{6vi{69f) = 2kVu{0^) - 4,^(Vi^Vi^) - eijknjVi{lk{S9f) 



(161) 
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Now, owing to the anisotropy caused by rapid rotation, we may write {Sv{59)'^) = {Sv±(S9)'^)l±/l± + {5vz{56Y)lz/lz 
and as {{56)^) is proportional to terms quadratic in and Iz, in the hmit k ^ and smah scales, one can easily 
reach at the following relation: 

{Sv^{5ef)\i^=o = Q (162) 

This relation predicts that in the presence of rapid rotation, and hence anisotropy, on the small line segment parallel 
to axis of rotation the correlation in the L.H.S. of (162) vanishes. This may be readily used in numerics to check if 
the two-dimensionalisation has been achieved and hence may be treated as a signature of the effect. 

V. GOY TURBULENCE 

In this section, we shall use GOY shell modcl[51, 52] (modified appropriately) to investigate the behaviour of the 
structure functions and, thus, the signatures of two-dimensionalisation effect. One may ask immediately why one 
needs another shell model though Hattori et. a?. [14] have already proposed a shell model — an improved version 
of shell model by L'vov couples of years back. To answer this question, let us collect the main results of that 
model: i) the exponent of the energy spectrum in the inertial range changes from —5/3 to —2, ii) no inverse cascade 
is detected with the increase in rotation rate, and iii) the PDF's of the longitudinal velocity difference doesn't match 
with the experiments. Well, this field of studying the two-dimensionalisation effect is growing rapidly. It has been 
confirmed well beyond doubt that the exponent overshoots the value —2 quite comfortably. One may refer to the 
experiments by Morize et. a/. [9]. The model justifies its results by invoking weak- wave-turbulence-theory in which 
inverse cascade is not really shown. This theory is a highly successful theory but one must be open-minded while 
dealing with problems as complex as turbulence and therefore, should take the experimental results at their face value. 
That some experiments and numerics do show inverse cascade with increase in the rotation rate should motivate one 
to construct shell models that can mimic this effect. As mentioned above, Hattori et. al.'s model finds PDF which 
mismatches with experiments and also, it requires a fluctuating part in the rotation rate to arrive at various results 
while in experiments and numerics there's no such part. This again should make it clear that why at all we need 
another model. Moreover, the numerical experiments done here are for unforced turbulence whereas Hattori et. aUs 
model dealt with forced turbulence. Hence, with due respect to the Hattori et. al.'s work, in this section we have 
tried to look at other possible shell model that can mimic the signatures of the two-dimensionalisation effect more 
closely. 

A. The model 

We have adopted the following strategy[53, 54] for the numerical experiments. A specific form of GOY shell model 
for non-rotating decaying 3D turbulence is: 

(163) 

This may be thought as a time evolution equation for complex scalar shell velocities u„(A:„) that depends on fc„ 
the scalar wavevcctors labeling a logarithmic discretised Fourier space {kn = ko2"). We choose: ko = 1/6,;^ = 
10~^ and n = 1 to 22. The initial condition imposed is: u„ = fc^/^e*^" for n = 1,2 and u„ = k^/'^e~''^e'^^" for 
n = 3 to 22 where 9n G [0, 27r] is a random phase angle. The boundary conditions are: w„ = forn < 1 and n > 22. In 
the inviscid limit {v — » 0), equation (163) owns two conserved quantities viz., J2n (energy) and X]ra(~l)"^n|^raP 
(helicity). If the fluid is rotating then one may modify equation (163) by adding a term i?„ = —i[u) + {—l)'^h] Un in 
the R.H.S. uj and h are real numbers. It may be noted that this term, as is customary of Coriolis force, wouldn't 
add up to the energy. The (— l)"/i term part in _R„ has been introduccd[15] to have non-zero mean level of helicity 
that otherwise has a stochastic temporal behaviour and zero mean level. Therefore, the appropriate shell model for 
rotating 3D turbulent fluid is: 

- i[uj + {-iyh]un (164) 

We fix /i = 0.1 in our numerical experiments and test for w = 0.01,0.1, 1.0 and 10.0. We shall henceforth refer w as 

rotation strength. Numerical results are obtained using 500 independent initial conditions and 40 different statistically 
independent runs. Inertial range has been taken as n = 4 to 15 — the range we are interested in. We have, by the 
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by, adopted slaved second order Adani-Baslifortli scheme[55] to integrate equations (163) and (164). 
The pth order equal time structure function for the model has been defined as: 

~ k-^" (165) 

Such has been done to avoid possible existence of period three oscillations[56]. The energy spectrum has been defined 
as: E{kn) = T,p{kn)/kn ~ k~"^. The mean rate of dissipation of energy is, of course, e = (X)„ J^fc^lwnP) and flux 
through nth shell is calculated using the relation: 

(166) 
(167) 

For studying relative structure function scaling, the ESS scaling exponents[57] are taken as = Cp/Ca- 't^j Cp and (* 
have all been calculated for inertial ranges only. 

B. The results 

The results are illuminating. (Detailed discussion and figures arc reported elsewhere [26]). One of them shows that as 
the rotation strength increases, the energy spectrum becomes steeper and the slope monotonically rushes from a value 
~ —5/3 to a value of ~ —7/3; hence validating one of the two-dimensionalisation effect's signatures. As we investigate 
into the direction of the flux in the inertial range regime, we can find that with the increase in rotation strength, first 
the forward cascade rate starts decreasing and then instances appear when at certain shells the flux direction reverses. 
Again, the number of such shells increase as the rotation strength is enhanced; clearly suggesting that depletion in 
the rate of forward energy cascade. Thus, yet another signature of two dimensionalisation has been upheld by the 
shell model. At this point, it must be appreciated how important the inclusion of term —i{—l)"h in equation (164) 
is in getting the effect of depletion in the rate of forward cascade. By setting mean level of helicity above zero, it is 
this very term that — in accordance with the arguments [24] that it is the helicity that is causing this signature of two 
dimensionalisation effect to show up — has empowered the model with the capacity to mimic the effect. Attempts 
to get this very effect by setting h — have failed miserably in our numerical experiments. The study of ESS in the 
shell model has been equally revealing. It has been noted that the increase in the rotation strength is accompanied by 
a departure from the usual She-Leveque scaling [58]. But, the fact that at higher p, (p seemingly becomes parallel to 
p/2 vs. p, is worth paying attention: This is in accordance with the direct numerical simulation (DNS) results[16] and 
experimental results[7]. However, most interesting observation would be that, within the statistical error, (* obtained 
for the rotating system via ESS coincides with that for the non-rotating ones. Probably, this extends the ESS for 3D 
fluids even further by implying that rotation keeps ESS scaling intact, even though usual Cp changes owing to rotation. 
Of course, only experiments and DNS can judge if this really is true for real fluid turbulence: GOY shell, after all, 
is just a model that remarkably reproduces many characteristic features of turbulence by only using a fraction of the 
computation power needed by DNS. In this context, one might be well aware that some modified versions of GOY 
model invented to model the distinguishing features of 2D turbulence have been shown to be rather useless[59]. One, 
thus, always has to be careful while dealing with simplified models of turbulence. 

VI. DISCUSSIONS AND CONCLUSIONS 

The Kolmogorov-Landau approach has been invoked in 2D homogeneous isotropic unforced fluid turbulence to arrive 
at the various correlation functions earlier obtained using different methods. Also, some experimentally verifiable 
correlation functions in the dissipation range have been derived. The results derived here are 'exact' (though not 
rigorous) something which is a far cry in the literature on turbulence. However, we have been careless enough to 
assume the existence of rj when ^ ^ 0. If 77 doesn't exist, the one-eighth law is in jeopardy. It is really unfortunate 
for the law that it has been rigourously proved [60, 61] that enstrophy dissipation is not possible for any 2D Euler 
solutions with finite enstrophy. Thus, i] may exist in the inviscid limit only when one takes rather ill-defined initial 
conditions for which the total initial enstrophy is infinite. In view of this, one must not be surprised at all if numerics 
and experiments fail to uphold the one-eighth law in many a situation. This very law of 2D turbulence, therefore, 
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doesn't enjoy the same classic status as the Kolmogorov law of 3D turbulence. 

Studies with structure functions of QG turbulence have again showcased how handy and useful the Kolmogorov- 
Landau approach can prove to be. The results (84) and (104) naturally agree with what has been arrived at by 
Lindborg[37] earlier. Within the domain of the approximations made these results are exact, something worth getting 
as the literature of turbulence is comparatively barren as far as exact relations are concerned. However, the hypothesis 
of the equipartition of energy used in equation (103) is as questionable as the assumption of isotropy in the sense 
of Charney. This hypothesis needs to be put on more firm basis. Again, the existence of Sq, like r], is questionable 
when u ^ 0. However, studies in 2D and QG turbulences in the perspective of the velocity structure functions 
hint that finding the velocity structure functions in the rotating flows with a view to unfolding the mysteries of 
two-dimensionalisation effect might not be just a wild goose chase. 

We emphasis on the fact that the form of two point third order structure function in a slowly rotating homogeneous 

3D turbulence can strongly indicate the initiation of the effect of two-dimensionalisation of 3D turbulence. It barely 
needs to be mentioned that the relations are quite interesting and pertinent (at least within the approximations made 
in the calculations) - something which, as is being said again and again, is worth getting in the literature of turbulence 
since exact relations are very few therein. So any theory developed in the limit of i?o and i?e — > oo, must satisfy 
the relation (149) derived in this article in the limit of low fl or explicitly violate the assumptions made to arrive at 
the result; in this sense the relation may prove to be of high importance. Moreover, true reason behind the so called 
two-dimensionalisation of turbulence has been figured out which accounts for the correct energy cascade direction 
and the correct energy spectrum found in the experiments and simulations. To settle the problem more neatly, it has 
been proposed that the study of passive scalars in rotating turbulence may prove to be of benefit. 
Shell models have been successfully used to study statistical properties of turbulence of many authors (see ref-([62]) 
and ref-([63]) for details). Most of them dealt with the homogeneous isotropic turbulence. Hattori et. al. gave a 
shell model for rotating turbulence. Their result could be improved in the face of the rapid progress in the field. We, 
thus, have used a modified version of GOY shell model to study the two-dimensionalisation effect. Some results of 
the model are, no doubt, consistent with experiments and DNS. However, one can always question the effectiveness 
of the signatures discussed herein because i) a scaling law for a single-component spectrum, though heavily used in 
literature, has poor meaning in the strongly anisotropic configuration relevant to pass from 3D-2D; different power 
laws can be found in terms of kz, k±_ and k in contrast to the 3D isotropic case, and ii) the inertial wave-turbulence 
theory is not consistent with an inverse cascade. Actually in weak- wave turbulence, getting rid provisionally of helicity 
and polarisation spectra, a two-component energy spectrum e(fc,cos0) with cos^ = k^/ \J (k1 + fc^) is found to be 
useful; if E denotes the traditional spherically averaged spectrum, the anisotropic structure is one of the best ways to 
quantify all intermediate states from isotropic 3D (with e = iJ(fc)/(47rfc^)) to 2D state (with e = Eikx^jji^-Eks^bi^z)). 
Two-dimensional trend is therefore linked to a preferred concentration of spectral energy towards the transverse wave- 
plane kz = 0. This concentration, however, does not necessarily yield an inverse cascade. A reasonable suggestion, 
in the light of this discussion, would be that in the shell model for rotating turbulence k should be interpreted as 
fc_L. It may be concluded that this study has put the equation (164) as a very firm shell model for the rotating 3D 
turbulent flows; after all, it explains the observed signatures of the two-dimensionalisation effect so closely. Probably, 
this model discussed in this article and the model due to Hattori et. al. can together yield a much better shell model 
for the rotating turbulence. 

In the absence of solutions of the Navier-Stokes equations, these various results regarding two-dimensionalisation can 
only be checked by data from experiments, and this endeavour seems to be difficult at present. So in the closing, we 

hope that validity of the results derived will be checked both numerically and experimentally in near future. After 
all, the structure functions traditionally provide checks for any plausible theory of turbulence. 
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